A dynamical theory of homogeneous nucleation for colloids and macromolecules 
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Homogeneous nucleation is formulated within the context of fluctuating hydrodynamics. It is 
shown that for a colloidal system in the strong damping limit the most likely path for nucleation 
can be determined by gradient descent in density space governed by a nontrivial metric. This 
is illustrated by application to low-density /high-density liquid transition of globular proteins in 
solution where it is shown that nucleation process involves two stages: the formation of an extended 
region with enhanced density followed by the formation of a cluster within this region.. 



Introduction Homogeneous nucleation is a fundamen- 
tal physical process of importance in fields as diverse as 
chemistry, materials science, biology and cosmology. Our 
basic understanding of it goes back to GibbshJ. The 
physics is governed by the fact that the excess free en- 
ergy of a liquid cluster relative to the vapor has a negative 
contribution that scales as the volume and a positive con- 
tribution due to surface tension that scales as the surface 
area. In Classical Nucleation Theory (CNT) it is as- 
sumed that the cluster is spherical, that its interior is in 
the bulk liquid state and that the surface tension is the 
same as for the coexisting liquid and vapor so that the 
free energy of the cluster as a function of its radius can be 
calculated giving a quantitative picture of homogeneous 
nucleation[l[. 

This description has several shortcomings. The surface 
tension is generally not constant and, especially for small 
clusters, can depend strongly on the size of the cluster. 
There is also no reason that the density within the cluster 
should be constant or, even if it were, that it should be 
equal to the bulk density for all cluster sizes. A more de- 
tailed picture can be devloped using Density Functional 
Theory (DFT) from which provides models of the free 
energy as a functional of the density profile(see e.g. Ref. 
[3,i). One can describe the density according to some 
parametrization (such as a hyperbolic tangent) which will 
involve at least three parameters: the central density, the 
radius and the interfacial width, and proceed by choosing 
a reaction coordinate - such as the radius of the cluster 
- and minimizing the free energy while holding the reac- 
tion coordinate constant. This does indeed lead to finite 
interfacial widths and size-dependent central densities, as 
expected(for recent examples, see Refs.[3, [H). 

Despite being physically reasonable, there are signifi- 
cant conceptual problems with this approach such as the 
arbitrarity of the reaction coordinate. As described, the 
nucleation pathway will consist of a monotonically in- 
creasing radius with the other parameters determined by 
the minimization. However, one could just as well choose 
the number of molecules in the cluster as the reaction co- 
ordinate in which case - in principle - the radius need not 
increase monotonically along the pathway since the mass 
of the cluster can increase by increasing the width while. 



at the same time, decreasing the radius 0. Even the ex- 
cess mass is not a good reaction coordinate in general 
and imposing more complex constraints can, at least for 
some models, lead to spurious divergences 

The fundamental difficulty underlying these and other 
equilibrium, free-energy based approaches is that the 
physical description is incomplete since homogeneous nu- 
cleation is a fundamentally nonequilibrium, fluctuation- 
driven process. This raises several questions about the 
classical description such as whether the free energy plays 
such a central role and whether it is necessary that the 
growing cluster actually pass through the critical clus- 
ter. The solution is to develop a nonequilibrium, dy- 
namical description of homogeneous nucleation and this 
is the goal here. The following development is based on 
Brownian Dynamics wherein molecules move according 
to Newton's laws while being subject to a frictional force 
as well as fluctuating forces. This is a simple model for 
colloids and the important case of macromolecules in so- 
lution in which case the friction and the fluctuations come 
from the bath/solvent. 

Theory The system consists of a collection of 
molecules of unit mass with positions and momenta qi,Pi 
interacting via a potential V. Additionally, the particles 
interact with a bath/solvent of light particles and this is 
described via a frictional drag and a fluctuating force 



dV 

<it = Pt, Pt = -T, 7Pi-H£j(t) (1) 



where all components of the fluctuating force are Gaus- 
sian and independent. 
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Defining the local density and momentum density respec- 
tively as 
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one sees that these satisfy the exact equations 
dp{r;t) 



dt 
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Coarse graining in space and assuming local equilibrium 
leads to a mesoscopic description of fluctuations in terms 
of fluctuating hydrodynamics. Neglecting temperature 
fluctuations gives 
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where p (r) and j (r) arc the coarse-grained local den- 
sity and momentum density, F [p] is the coarse-grained 
free energy and 11 is the dissipative part of the stress 
tensor which has both a deterministic and a fluctuating 
contribution 0. The free energy term is a representation 
of the local pressure and has been discussed extensively 
in the DFT literature: its use here can be viewed as a 
local equilibrium approximation!^ [sl-fioj. The quantity 
^ (r; t) is the noise due to the Brownian dynamics and is 
white and delta-correlated in space and time. Note that 
this is just the natural generalization of Landau and Lif- 
shitz's fluctuating hydrodynamics taking account of the 
Brownian forces. Assuming that the velocity will always 
be small due to the damping, the convective term can 
be neglected so that the second equation becomes linear 
in the momentum density. Eliminating the momentum 
current then gives 
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In the following, the second-time derivative, the so-called 
incrtial term, will be neglected, as is usual in the strong- 
damping approximation. Then, when the density is low, 
in the ideal gas limit, the first term on the right becomes 
7~^fcsTV^p(r) so that D = ^~^kBT can be identified 
as the diffusion constant. 

The use of fluctuating hydrodynamics as basis for 
studying nucleation is similar to the approach devel- 
oped by Langer[llj. The primary difference here is that 
the emphasis is on understanding the time-evolution of 
the formation of the critical cluster whereas previous 



work focused on the nucleation rate. This development 
differs from more phenomenological approaches which 
arc couched entirely in terms of order parameters, such 
based on nonequilibrium thermodynamics [lij or phase 
field theory [I^, in that nonlinearities of the transport 
coefficients and colored noise occur naturally and play 
an important role. One of the goals below is to relate the 
hydrodynamic description to one involving order param- 
eters. 

In order to characterize the generic properties of the 
process of nucleation, we focus here on the most likely 
path (MLP) where a "path" is understood as a func- 
tion p(r; t) connecting the initial state of pure metastable 
phase and the final state of pure stable phase. When the 
noise amplitude is small (as in the strong damping limit), 
most systems should go through a nucleation pathway 
close to this generic result. In general, determining the 
MLP is complex. However, without the inertial term, 
Eq.© is a gradient-driven, diffusive dynamics which 
obeys a fluctuation-dissipation relation. By a straightfor- 
ward generalization of [l4| , it can be shown that for this 
type of dynamics the MLP connecting metastable states 
does indeed pass through the critical point and that it 
coincides with either the forward-time or backward-time 
deterministic trajectory in density space[l^. The MLP 
can therefore be determined by starting at a local mini- 
mum and moving along the deterministic path 
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where the sign is chosen according to the direction one 
wishes to movejl4 Il5|. 

Equation® is the primary theoretical result of this pa- 
per. It superficially resembles the usual Dynamic Density 
Functional Theory (DDFT) equationdHi,!!^] but is in 
fact considerably more general. It says that the most 
likely path can be determined by following the DDFT 
dynamics when that dynamics does indeed connect the 
desired states such as in passing from a high-energy to 
a low-energy state with no barrier separating them. (An 
example of this would be spinodal decomposition[9|.) 
However, DDFT cannot describe the crossing of a free 
energy barrier as it specifically pertains to the ensemble- 
averaged density. In essence, it is the result of averaging 
Eq. ([7]) (without the inertial term) over the noise. In con- 
trast, Eq. ^ also describes the MLP when this means 
going uphill against the free energy gradient. In that 
case, it says that the MLP can be obtained by revers- 
ing the sign of the gradient or , equivalently, by following 
the time-reversed dynamics [IJ, |l5[. It can therefore be 
viewed as an extension of DDFT to barrier-crossing prob- 
lems, given the various assumptions set out above. This 
simple result is strongly dependent on the existence of 
the fluctuation-dissipation relation in Eq.Q and on the 
assumption of weak noise (compared to the thermody- 
namic driving force) . It will not be exact if either of these 



conditions are violated and, in particular, the much more 
complicated strong-noise result will be discussed in at a 
later time. Finally, it is important to note that Eq.® 
is simply a mathematical means of identifying the MLP 
and that it does not imply in any sense that the actual 
(strongly-dissipitive) dynamics is time-reversal invariant. 

In principle, Eq.® could be integrated directly to de- 
termine the MLP or some other technique, such as the 
string method, used to determine the path. However, 
the goal here is to generalize previous descriptions of nu- 
cleation which are based on a set of order parameters 
characterizing the system. In CNT, the cluster is as- 
sumed to be spherical and only parameter is the size of 
the cluster: more generally, a minimal set would include 
some measure of the density inside the cluster and the 
width of the interface as well. In the present formal- 
ism, the order parameters must somehow be related to 
the spatial density since it is the fundamental quantity 
describing the evolution. We therefore imagine that the 
density profile is approximated by some test function of 
the form p{r,t) = /(r;x(t)) where x (i) stands for the 
set of order parameters. It is possible to give an exact 
equation for the evolution of the parameters based on an 
analysis of the MLP but here a more heuristic method is 
used. First, Eq.([8|) is integrated over a spherical volume 
of radius r giving 
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where m{r;x{t)) is the mass in side the spherical shell 
of radius r and the notation indicates a surface integral 
over that shell. Then, spherical symmetry is assumed and 
Eq.® is multiplied by (r) ^^^^ and integrated 

over r to get 
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where the metric is 
gab = 47r 
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and where O = F—p,N is the grand potential which arises 
due to an integration by parts. This becomes exact if the 
parametrization is complete in the sense that / (?-; x (t)) 
is able to represent any well-behaved function arbitrarily 
closely(e.g. an expansion in a complete set of basis func- 
tions). The exact minimization of the action for the case 
of a finite number of order parameters and its relation to 
this approximation will be discussed elsewhere. 

Application to low density /high density liquid transi- 
tion in globular proteins Many proteins in solution ex- 
hibit a phase transition between a low density gas-like 
phase and a high density liquid-like phase. This behavior 
can be modeled using an effective pair-potential in which 
case it becomes analogous to the vapor/liquid transition 




FIG. 1. The excess particle number, excess free energy and 
central density relative to their values in the critical cluster 
as functions of the natural reaction coordinate, s, where s=0 
corresponds to the initial vapor phase and the final liquid 
phase occurs for s = oo. The vertical line marks the transition 
between the two growth regimes (see text). 



in simple fluids. Calculations were performed for the ten 
Wolde-Frcnkel[l6l| model potential for globular proteins 
having hard-core diameter a and energy scale e using the 
squared-gradient free energy model. 



/(p(r)) + \k{Vp{v)Y ) dv (12) 



F[p] 



where f{p) is the bulk free energy per unit volume, cal- 
culated using thermodynamic perturbation theory, and 
the coefficient K was calculated using a recently-derived 
approximation^. Equation ([S]) was integrated (assum- 
ing spherical symmetry) by discretizing the right hand 
side using Ea.([TO| and the method of piece- wise linear 
profiles 13 (equivalent to a variable grid method). At a 
temperature of fc^T = 0.375e and with a pressure su- 
persaturation of 1.159, the exact excess energy barrier 
was found to be A/3f2 = 75.8 with 1158 molecules in the 
cluster while the discretization with 19 parameters gives 
a value of 77.1 and 1175 molecules. The MLP was then 
determined by starting near the at the critical cluster and 
perturbing in the direction of the ne gati ve eigenvalue [l7|| 
and integrating Ea. (|10p numerically [iSj. In tracing the 
backwards part of the path, the calculations were termi- 
nated when /3Af2 = IksT since the weak noise approxi- 
mation is not applicable for lower energies [lij. 

Figure [T] shows the evolution of the excess number of 
particles, excess free energy and of the central density. 
The independent variable is the natural reaction coor- 
dinate which is distance along the nucleation pathway 
as calculated using the metric, Ea.([TT|). When the clus- 
ter is large, the path is similar to that which would be 
obtained using typical heuristic methods. However, for 
smaller droplets, the results are quite nonclassical. Fig- 
urc[2]gives the spatial size of the droplet according to two 
different measures: the equimolar radius as calculated 
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FIG. 2. The same as Fig. [T] but showing the equimolar radius, 
Req, and the total spatial extent, Rtotai, along the path. The 
vertical line marks the same point as in Fig[T] 



based on the central density and the (model-dependent) 
total spatial extent of the droplet (in this model, the 
droplet always has a well-defined finite support). Com- 
bining the information in these two figures, it is seen 
that the MLP begins with a spatially-extended distur- 
bance having very low density but a fixed excess number 
of molecules (in the present case, about 450). 

The interpretation of these results is not as different 
from the usual picture of nucleation as they might at 
first appear. At short times, during which the equimo- 
lar radius is nearly constant, a small increase in density 
forms over a spatially extended volume. From Fig. 1, 
it is apparant that, despite its spatial extent, the excess 
energy of this density fluctuation is quite small so that 
its formation is not improbable. The second part of the 
process is the formation of a cluster within this region of 
enhanced density. That nucleation would preferentially 
take place in a region of enhanced density, which there- 
fore already contains some of the excess mass needed to 
form a cluster, is quite reasonable. Indeed, having the 
necessary excess mass present locally allows the cluster 
to form more rapidly than if matter had to diffuse in 
from the surrounding bulk. What does appear strange 
is the directed nature of process with the density fluc- 
tuation appearing to contract to form the nucleus. This 
is partly a result of insisting on spherical symmetry and 
partly due to a well-known property of the MLP wherein 
it typically involves a system crossing a barrier in the 
shortest number of steps possible with no back-tracking 
or variation[l^. The MLP is of course an abstraction: 
an actual realization of the processes will involve fluctu- 
ations around it and will not appear so deterministic ( see, 
e.g., Ref [Tgji ) . Finally, because mass is conserved, any dy- 
namical process of cluster formation is going to give the 
appearance of drawing in mass from the boundaries of 
the system. 



Conclusions A description of nucleation applicable to 
colloids and macromolecules in solution based on fluctu- 
ating hydrodynamics has been developed. It was shown 
that under assumptions of strong dissipation and weak 
noise the most likely path could be determined by gra- 
dient descent on the free energy surface and that it nec- 
essarily passes through the critical point, thus provid- 
ing justification for more heuristic methods based solely 
on free energy considerations @. It is also interesting 
to note that Eq. pH)) can be seen to justify more phe- 
nomenological treatments of nucleation in which a set 
of order parameters is assumed to evolve stochastically 
as dx/dt = i-§f- + ? with a fluctuation-dissipation rela- 
tion determining the amplitude of the noise. The same 
approach can be applied to other nucleation phenomena 
such as heterogeneous nucleation, nucleation in confined 
systems and even, conceivably, to transitions in granular 
fluids. 
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